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1. Introduction 

Quantum Chromodynamics (QCD), the theory of strong interactions is a non-Abelian quantum 
field theory characterized by local SU(A'^) gauge invariance where N = 3 denotes the number of 
"colours". The adjoint gauge bosons (gluons) couple nf "flavours" of fermionic matter fields in the 
fundamental representation (quarks). 

The standard non-perturbative definition of QCD is based on lattice regularization [1], which 
makes the theory mathematically well-defined and amenable to analytical as well as numerical 
studies. Theoretical progress, algorithmic innovation and ever more powerful computers have al- 
lowed some teams to derive non-perturbative low energy properties like the spectrum of QCD with 
nf = 2, rif = 2 + 1 and = 2 + 1 + 1 sea quarks at realistic values of the physical parameters, see, 
e.g., ref. [2] for a recent review. 

A different non-perturbative approach to QCD is based on an expansion in powers of of 
the inverse number of colour charges at fixed rif [3]. When N is taken to infinity, and the gauge 
coupling g is sent to zero, keeping the product g^N as well as nf fixed (the 't Hooft limit), the 
theory reveals striking mathematical simplifications, see refs. [4] for recent reviews. For instance, 
all amplitudes of physical processes are determined by a particular subset of Feynman diagrams 
(planar diagrams), the low-energy spectrum consists of stable meson and glueball states and the 
scattering matrix becomes trivial. One may study the physical N = 3 case, expanding around the 
large-N limit in terms of 1 /N. Interestingly, the non-flavour- singlet spectra of QCD with sea quarks 
and quenched QCD agree within 10 % [2], which may indicate both nf/N and 1 /N^ corrections to 
be small. 

The large-A/^ limit also plays a vital role in the chiral effective theory approach where the 
A'^-dependence of low-energy constants is known [5] and, within this framework, in studies of 
properties of unstable resonances, see, e.g., refs. [6, 7]. (Un)fortunately, even in the large-A^ Umit 
QCD is far from trivial. 

Another non-perturbative approach — that, unlike lattice regularization, does not break the Eu- 
clidean spacetime symmetry — to low-energy properties of strongly coupled non-Abelian gauge 
theories is based on the conjectured correspondence between (supersymmetric) large -A'^ gauge and 
string theories in the classical gravity limit in an anti-de-Sitter spacetime (AdS/CFT correspon- 
dence) [8]. During the last decade, many studies have used techniques based on this coiTcspon- 
dence to construct models which reproduce the main features of the meson spectrum of QCD [9]. 

Recently, the dependence of various quantities on A'^ was studied in lattice simulations. For 
instance, pseudoscalar and vector meson masses (among other observables) were determined in 
refs. [10-13]. Here we improve upon and extend these studies, reducing the quark masses and 
increasing N, the statistics, the number of states studied and the volume. 

In view of the above discussion, it is important to determine the meson spectrum of large- 
QCD to constrain effective field theory parameters and to enable comparison with AdS/CFT 
and AdS/QCD predictions. We perform our simulations in the quenched approximation to QCD, 
neglecting sea quark loops. Therefore, we only encounter 1/A'^^ corrections to the large-A'^ limit, 
rather than nf/N coiTcctions. This allows for a more constrained N ^ °° extrapolation, at the 
same time reducing the computational effort. We remark, however, that the naive cost of including 
sea quarks into the update only scales like A'^^ while the pure gauge operations scale like A'^^ : in 
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the large-A'^ limit not only the quenched theory becomes unitary and identical to the un-quenched 
theory but so also does the computational effort, which is quite substantial &tN = 17. 

Finally, we aim at clarifying a discrepancy between the results of refs. [10- 12], which at large 
favour a value of the vector meson mass close to that of real-world QCD, and those obtained in 
ref. [13], reporting a value approximately twice as large. 

2. Simulation details 

Our simulation strategy is to tune the lattice couplings, keeping the square root of the string 
tension y/aa « 0.2093 in lattice units a fixed. We employ the Wilson gauge and fermionic action 
and aie not yet in the position to perform a continuum limit extrapolation. Lattice artefact terms will 
have the same functional large-A'^ scaling as the dominant continuum limit terms and hence basically 
should not affect the size of corrections. Our experience with the present action [10, 11] leads 
us to expect systematic errors on mass ratios of about 5 % at the present lattice spacing. 

We define the decay constant of a meson X in the large A'^ limit as 




(2.1) 



where Fx = fx/ \/2- We distinguish between « 92 MeV at physical quark masses and F = 
Fniniq = 0). 

Using the ad hoc value a = 1 GeV/fm, our lattice spacing corresponds to a ^ 0.093 fm or 
a^' S3 2.1 GeV. Strictly speaking, we can only predict ratios of dimensionful quantities. In the real 
world where experiments are performed, Uf >0,N = 3 and even the string tension is not well 
defined. This means that any absolute scale setting in physical units will be arbitrary and is just 
meant as a rough guide. Nevertheless, we notice that other ways of setting the scale give similar 
results, indicating that the N = °° world is not far removed from N = 3 QCD with sea quarks. For in- 
stance, in the chiral limit we find F°° = 0.22(2) = 96(9) MeV, in qualitative agreement with the 
real QCD value [14] F = 85.9(1.2) MeV. Moreover, we obtain trip = 1.638(7)^^ = 728(3) MeV 
at large A'^, quite close to the experimental p -meson mass of 775 MeV. This is remarkable, in par- 
ticular since this resonance has a decay width of almost 150 MeV. 

The string tension was computed in ref. [15] for A/^ = 2,3,4,6 and 8 and in ref. [16] for A^ = 5 
and 7. For A'^ = 17 no value is known and our lattice volume (see below) is too small for a reliable 
determination from torelon correlators. Therefore, we estimate 

A w a^^exp 

and extrapolate the A-parameter obtained at A'^ < 8 for a^/a = 0.2093 as a polynomial in 1 /N^ to 
A^ = 17. The lattice coupling is defined as jS = 2N^/X = N/{2na), where A = A^^^ is the 't Hooft 
parameter in the lattice scheme. We obtain the central value /3i7 = 208.45^29 from a 1/A'^^ fit to 
the N >6 data, with systematics estimated by varying the fit range and allowing for a quartic term. 
Our j8 -values and simulated volumes are summarized in table 1. Note that the lattice 't Hooft 
couplings deviate by 1/A'^^ terms from a constant along our trajectory of fixed aa^. We remark 
that incidentally ref. [13] simulated SU(17) at j8 = 208.08 which is almost identical to our value 
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N 


N^, X N, 




A 


lO^K- 


'^conf 


2 


16^ X 32 
24^ X 48 
32^ X 64 


2.4645 


3.246 


14581, 14827, 15008, 15096 
14581, 14827, 15008, 15096, 15195.9 ,15249.6 
14581, 14827, 15008, 15096, 15195.9 ,15249.6 


400 
200 
100 


3 


16^ X 32 
24^ X 48 
32^ X 64 


6.0175 


2.991 


15002, 15220, 15380, 15458 
15002, 15220, 15380, 15458, 15563.8, 15613 
15002, 15220, 15380, 15458, 15563.8, 15613 


200 
200 
100 


4 


16^ X 32 
24^ X 48 


11.028 


2.902 


15184, 15400, 15559, 15635 
15184, 15400, 15559, 15635, 15717.3, 15764 


200 
200 


5 


16^ X 32 
24^ X 48 


17.535 


2.851 


15205, 15426, 15592, 15658 
15205, 15426, 15592, 15658, 15754.8, 15835.5 


200 
200 


6 


16^ X 32 
24^ X 48 


25.452 


2.829 


15264, 15479, 15636, 15712 
15264, 15479, 15636, 15712, 15805.1, 15884.5 


200 
200 


7 


16^ X 32 
24^ X 48 


34.8343 


2.813 


15281.6, 15496.7, 15654.7, 15733.9 
15281.6, 15496.7, 15654.7, 15733.9, 15827.3, 15906.2 


200 
200 


17 


12^ X 24 


208.45 


2.773 


15298, 15521, 15684, 15755, 15853.1, 15931 


80 



Table 1: Simulation parameters: ?iconf is the number of configurations analysed for eacli set of 
parameters. All configurations were separated by 200 combined heatbath and overrelaxation Monte 
Carlo sweeps and found to be effectively statistically independent. 



j3i7 = 208.45, ruling out the hypothesis that our N <\1 simulation points are not in the "continuum 
phase of the large-A^ theory" [13]. In fact at our lattice spacing we will find the 7 > N >3 spectra 
to be in almost perfect agreement with the N = 17 results. 

Different lattice sizes are investigated to exclude finite volume effects spoiling the large-A'^ 
extrapolation, in particular at the lighter quark masses. As expected, these become iiTclevant at 
large A'^, thereby justifying the use of a relatively small lattice atN = 17. 

The so-called hopping parameter k is related to the lattice quark mass via 

Kc denotes the critical value, corresponding to a massless quark. The additive constant is given by 
K^^ = 8 + ^( A ) and we will determine this non-perturbatively. The fc- values shown in table 1 were 
selected to keep one set of pion masses approximately constant across the different SU(A'^) theories. 
We vary the "pion" mass down to niji/ y/a 0.5 for groups with N >5, and to m^j/ ^/a « 0.75 for 
N < 5. We also simulated a smaller quark mass for SU(A'^ < 5) but found significant numbers of 
"exceptional configurations" [17] (up to 15 % of the total); we leave these data out of this work. 
For N = 5,a.t the lowest quark mass, only two exceptional configurations were encountered that we 
removed from the analysis. 

Our code is based on the Chroma suite [18], which we have adapted to work for generic A'^ 
values. We compute correlation matrices between differently smeared interpolators, allowing us 
not only to extract the ground states but also giving us access to excitations in many channels. 
Details of the analysis can be found in ref. [19]. 
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3. Results 



We employ the lattice quark mass mpcAC> defined through the axial Ward identity, as our ref- 
erence mass, avoiding the additive renormalization (ff^r^) of the quark mass rtiq defined in eq. (2.3). 
These two masses are related by a combination of renormalization constants, 

Zp , , 1 / 1 1 \ 
a'WpCAC = 7;r^ (1 +o«»ipCAcH )^\ , (3.1) 

where the bam-teim parameterizes the leading lattice correction. All three fit parameters Kc, b 
and Zp/{ZaZs) are well described by constants plus l/A/^^ -corrections. We find the latter ratio 
of renormalization constants to vary between 0.68 (A'^ = 2) and 0.83 (A'^ = 17), with the SU(3)- 
value 0.75, which is consistent with the non-perturbative result 0.81(7) [20] obtained at j3 = 6.0 
— close to our value j8 = 6.0175. Motivated by the weak A'^-dependence of this result — which 
is also supported by perturbation theory — we will use the non-perturbative SU(3)-values for 
the renormalization factors of quark bilinears for all gauge groups, allowing for a 8 % systematic 
uncertainty, due to this approximation. These factors are needed to determine the decay constants 
below. 

Next, we determine the dependence of the pseudoscalar mass mn on the quark mass mpcAc: 



m 



A quenched chiral log is expected at small A/^- values, parameterized by 8. We also include a sub- 
leading term to prevent interference between the larger mass data and the chiral log. The fits for 
the different A'^ are depicted in figure 1. 5 is expected to be suppressed by a factor [5] 1 /N, with 
1/A'^^ corrections. We find the parameter values 

A i9 9Vnim 9-7(1.6) 6.3(2.2) 0.021(19) 1.12(21) 

A = 12.23(0.10) , B= 1.74(0.13) + — p—, d = + ^^3 — '(3-3) 

see the right panel of figure 1. 5 > 0.05 at A'^ = 3 is coherent with expectations but the 1 /N coeffi- 
cient is statistically compatible with zero. In principle, the determination of 5 may be obscured by 
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Figure 2: The spectrum, extrapolated to the chiral limit. The error bands correspond to N = °°. 

the possibility of a non-zero pion mass at mpcAC = 0, at finite lattice spacings. However, the linear 
disappearance of m| oc mpcAC for A'^ > 5 indicates chiral symmetry to be broken only mildly. 

We choose to parameterize the quark mass dependence of all our results through mpcAC 
which can be determined more precisely and reliably than the pseudoscalar mass niji. The relation 
m^(mpcAc) above enables the translation of a functional dependence on mpcAC into a dependence 
on nijc- Nevertheless, we quote the N ^ °° result 

1+0.375(64) f-^^ +•••, (3.4) 



mp{mjt) 



mp(0) ■ \mp{0)_ 

to allow for a direct comparison with the prediction, e.g., of ref. [21]. 

We display the chirally extrapolated spectrum in figure 2 for the different SU(A'^) groups. The 
— oo values are shown as horizontal error bands. Note that the N = 17 values are perfectly 
consistent with the N < 17 results, ruling out the twice as large p-meson mass obtained at almost 
the same lattice coupling in ref. [13]. This may be an artefact of the method used in that reference, 
e.g., due to excited state pollutions. 

The figure does not only give the masses of the excited and ground state mesons but also 
illustrates the decay constants /;,: and fp to show the expected scaling behaviour oc -y/^V, with small 
corrections. The analysis of the scalar (ao) correlation function at low A'^ is complicated by the 
presence of ghost states, due to the unitarity violation of the quenched model. Subtracting these 
contributions results in large errors. 

Of phenomenological interest is not only the spectrum at nig = but are also the spectra at 
mq = Mud and at iriq = where m,„/ and denote the physical (isospin-averaged) light quark and 
strange quark masses, respectively. We fix the former, imposing the value [14] 



1.073(15) 



(3.5) 



6 



The meson spectrum in large-N QCD 



Gunnar S. Bali 









Particle 




mq = 


mq = Mud 


Mq = 


Mq = 


mq = Mud 


mq = nis 


TZ 


+ 





0.417(100) 


1.62(10) 





1.92(46) 


7.46(48) 


P 


1 


1.5382(65) 


1.6382(66) 


1.9130(79) 


7.08(10) 


7.54(11) 


8.80(13) 


ao 


0++ 


2.401(31) 


2.493(31) 


2.755(32) 


11.04(21) 


11.47(22) 


12.67(23) 


ai 


1++ 


2.860(21) 


2.938(21) 


3.158(22) 


13.16(21) 


13.51(21) 


14.53(23) 


bi 


i+- 


2.901(23) 


2.978(23) 


3.197(23) 


13.35(21) 


13.70(22) 


14.71(23) 


TZ* 


0-+ 


3.392(57) 


3.462(57) 


3.659(58) 


15.61(34) 


15.93(35) 


16.83(36) 


P* 


1 


3.696(54) 


3.756(54) 


3.928(54) 


17.00(34) 


17.28(35) 


18.07(36) 


«0 


0++ 


4.356(65) 


4.420(65) 


4.603(66) 


20.04(41) 


20.33(41) 


21.18(42) 


a\ 


1++ 


4.587(75) 


4.646(75) 


4.816(77) 


21.10(46) 


21.38(46) 


22.15(47) 


b\ 


1+- 


4.609(99) 


4.673(99) 


4.85(10) 


21.20(54) 


21.50(55) 


22.33(56) 


Jk 




0.3074(43) 


0.3271(44) 


0.3784(56) 




1.505(29) 


1.741(36) 


POO 

Jp 




0.5721(49) 


0.5855(50) 


0.6196(64) 


2.632(43) 


2.694(44) 


2.850(50) 



Table 2: The N = °° meson spectrum and decay constants in units of the square root of the string tension a 
and in units of the (normalized) chiral pion decay constant Foe ~ F \/3/N for three different values of the 
quark mass. A systematic error of 5 % needs to be added, due to the missing continuum limit extrapolation. 
Because of the non-perturbative N — 3 rather than N = °° renormalization, an extra 8 % error should be 
added to the last three columns and to the last two rows of the table, with the exception of the /^/F"-ratios 
where this factor cancels. 

at A'^ = 3, keeping / ^/a constant for A/^ / 3. The renormalization constant and A/^-dependence 
cancel from the ratio. The strange quark mass is obtained by fixing the ratio of a (fictitious) strange- 
antistrange pion over the (1020) vector particle at N = 3 to the experimental value 

mnjms) _ 686.9 

mp(m,) ~ 1019.5 ' ^ ' ' 

where (m^± +m^o ~"J^±)^^^ ^ 686.9 MeV. We display the results at the different quark masses 
in table 2. The normalization in units of F°° should be particularly useful for chiral perturbation 
theory applications [6, 7]. 

4. Summary 

We have determined the decay constants as well as the ground and first excited state masses of 
mesons in the large-A'^ limit of QCD by lattice simulation of the A'^ = 2,3,4, 5,6,7 and 17 quenched 
theories. In almost all but the scalar channels 1 /N'^ corrections are found to be tiny for A'^ > 3. In 
a forthcoming publication [19] we will compare our findings to model expectations. We find the 
scalar to be by factors of about 1 .5 heavier than the vector particle at any quark mass smaller than 
the strange quark, which is of particular relevance to the phenomenology of scalar mesons. 
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